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حکیده 


تغییرات ارتفاعی رخداده در سطح زمین که غالبا ب‌صورت فرونشست و گاهی بالاآمدگی ظاهر می‌شود. 
به‌عنوان یک مخاطره خحاموش ولی جدی در محدوده شسهرها بهحساب می‌آید. که می‌تواند به‌مرورزمان 
صدمات جدی به تأسیسات شهری وارد سازد. به دلیل تغییرات رخداده در سطح آب‌های زیرزمینی شهر 
کرمان پدیده فرونشست و بالاآمدگی با شدت زیادی در این محدوده وجود دارد که می‌تواند خسارات 


جبران‌ناپذیری به شهر وارد نماید و هدف تحقیق پایش این تغییرات است. در این تحقیق سعی شده است 
که با استفاده از تکنیک تداخل‌سنجی تصاویر راداری میزان تغییرات ارتفاعی رخداده در محدوده شهر 
کرمان بررسی و نحوه تغییرات مکانی آن در طول ۱۶ سال اخیر مورد پایش قرار گیرد. در اين راستا از ٩‏ 
و ار رک وت اه ۱ طبار اور فان نتفای 


شد و با انجام تکنیک تداخل سنجی» چهار تداخل‌نگاشت از محدوده مورد مطالعه تهیه گردید. با مطالعه 


تداخل‌نگاشت‌هاء نرخ و دامنه فروند ت و بالاآمدگی استخراج گردید. بر این اساس حداکثر نرخ 
فرونشست و بالاآمدگی در چهار دوره زمانی مربوط به سال‌های ۱۳۸۲ - ۱۳۸۳ ۱۳۸۲-۱۳۸۹ ۱۳۹۱- 


سانتیمتر در سال بالاآمدگی بوده است. استخراج عرصه در معرض مخاطره نشان داد از مجموع مساحت 
محدوده در حدود ۶۳ درصد در پهنه‌های پرخطر تا نسبتاً پرخطر قرار دارد. شواهد میدانی نشان می‌دهد که 
علاوه بر محدوده‌های فرونشستی» در مناطق با نرخ تورم و بالاآمدگی زیاد نیز آثار و شواهد خسارات به 
ساختمان‌ها به‌وضوح دیده می‌شود. نقشه‌های جابجایی ایجاد شده نشان می‌دهد که سطح زمین در شهر 


۱ نویسنده مسئول: ٩۱۳۱۹۳۶۲۵7‏ ۵0۵۲ زحامطه ]۱۷ :1تظرظ 


۱۳3۹/4۹ جغرافیا و مخاطرات محبطی شماره سی‌ام 


کرمان از سال ۱۳۸۳ تاکنون دجار یک روند فرونشست فزاینده‌ای شده است به‌طوری که علاوه بر افزایش 
نرخ قرو تکششگان او فه‌هاعن رف ازشهن فر کی آن شاه است: 


۱. مقدمه 
یکی از مشکلات اساسی و روزافزون در اکثر جوامع بشری که بیشتر در اثر فعالیت‌های بشری رخ می‌دهد پدیده 
فروندٌ است. فرونشست به فرورفتگی سطح زمین اشاره دارد که بر اثر عوامل گوناگون طبیعی مانند انحلال» آب 


رشد روز افزون جمعیت منجر به استفاده بی‌رویه از آب در مصارف خانگی» صنعتی و کشاورزی شده و پیرو آن ایجاد 


اثرات نامطلوب کمی و کیفی در منابع آب بوده است. افزایش استفاده از آب‌های زیرزمینی به‌ویژه در حوضه‌هایی که با 


نهشته‌های آبرفتی» دریایی کم‌عمق یا دریاجه‌ای تحکیم نيافته انباشته گشته‌انده باعث ایجاد فرو می‌شود. در 


بسیاری از دشت‌های ایران بهره‌برداری بیش ازحد از آب زیرزمینی منجر به رخداد فروندٌ ِ در آنها شتاله اشسبت 
(یاراحمدی» ۱۳۹۵؛: 2016 .21 61 12]71). (شریفی کیاء ۱۳۹۱؛ صالحی و همکاران ۱۳۹۲؛ ,07067 0ص رد«مال0) 
1 2016 .21 02211870 2017 .21 61 «اش1). بدین‌جهت. احتمال دستیابی به شواهدی در ارتباط با وجود پدیده 


فرونشست در دشت‌ها بالاست. از آنجائی که اولین گام در مطالعه این پدیده. شناخت گستره فضائی آن و اندازه‌گیری 


هر جچه دقیق‌تر میزان فروندٌ 5 است؛ بنابراین با شناحت مشخصات مکانی و رفتار زمانی این پدیده» امکان ارائه و 


توسعه مدل منطقه‌ای این پدیده میسر شده و از این طریق می‌توان راهکارهایی اساسی و عملی در جهت جلوگیری از 


روند آتی و کاهش خحسارات مرتبط با آن اتخاذ کرد. امروزه. برای اندازه‌گیری فرونث ت زمین روش‌های مختلف 
ژودیتیکی و غیرژئودیتیکی وجود دارند. مشاهدات زمینی از طریق برداشت‌های توتال استیشن و لیزر اسکنر 


ایستگاه‌های 5 اندازه‌گیری‌های دقیق و پیوسته‌ای را در تقاط محدودی از منطقه فروندٌ ِ ارائه می‌دهند بنابراین 


هیچ کدام از این روش‌ها. توانائی تعیین وسعت و الگوی فضائی پدیده فروندٌ 9 را ندارند (2004 ,21 )6 اتمه 1 


02 هنهک ,2015 .له 64 20 از طرف دیگر تکرار انجام هرکدام از این روش‌ها به‌ویژه هنگامی‌که برای 
منطقه وسیعی مورد نظر باشد بسیار پرهزینه و زمان‌بر خواهد بود. وجود چنین محدودیت‌هایی همواره از چالش‌های 
اساسی محققان در ارتباط با اندازه‌گیری‌های دقیق و پایش فضایی تغییرات سطح زمین محسوب می‌شود. با ظهور 
سنجنده‌های دورسنجی مختلف و ماهواره‌های راداری تصویربرداری در دهه ۱۹۹۰ و توسعه سریع آن در سال‌های 
بعدازآن افقی تازه و بسیار امیدوارکننده‌ای پیش روی محفقان علوم زمین قرارگرفت (یاراحمدی» ۱۳۹۵). همگام با 


توسعه روش‌های مختلف پردازشی این قبیل داده‌های دورسنجی بکارگیری روش اینترفرومتری راداری " نیز در سال 
های اخیر به‌عنوان یکی از روش‌های غیرژئودیتک و ابزاری کارآمد در آشکارس ازی تغییرات سسطح زمین و پایش 
جابجائی‌های ناشی از پدیده‌های مختلفی نظیر فرونٌ ت» آتشفشان و زمین لغزش و زلزله و ... مورد توجه محققان 


علو م زمین قرار گرفته است (2004 .21 اه عصقط :2007 .21 اه طعه/۱۷]0 :2011 .له )ع وزانن؟ رتصتم‌عم؟ .له اه بان1 


4 :2016 ,2020 :2013 : 2017 .21 )6 هنعط مهرابی و پورحسروانی» ۱۳۹۷). برخی مزایای این روش که 
نسبت به نوع سنجنده متفاوت است می‌توان به پوشش مکانی وسیع و پیوسته امکان محاسبه جایجایی‌ها با دقت 
سانتیمتر و پائیش قابلیت تصویربرداری در هر شرایط آب و هوایی و هر ساعت از طول شبانه روز تکرار 
تصویربرداری از یک منطقه حاص در هر ۳۵ یا ۱۱ روز و سهولت دسترسی به داده‌های آرشیوی 5 آنلاین سنجنده‌های 
راداری و در نتیجه نیاز کمتر به انجام مجدد عملیات میدانی و صرف هزینه بسیار پائین آن در مقایسه با روش پیشین 
ذکر کرد (یاراحمدی, ۱۳۹۵). در دهه‌های اخیر. افزايش جمعیت جهان به خصوص در مناطق شهری به‌عنوان یک 
پدیده مهم پیچیدگی‌ها و مشکلات زیادی را در زمینه‌های مختلف ایجاد کرده است. در این میان با برداشت بی‌رویه از 
آب‌های زیرزمینی خطر فرونششست زمین در محدوده‌های شهری به‌عنوان یکی از مهمترین موضوعات در حوزه‌های 
برنامه‌ریزی و طراحی شهری به جهت مدیریت بحران و حطرپذیری. کاهش ریسک و آسیبها وافزایش ایمنی می‌باشد. 
مناطقی شهری به دلیل تراکم جمعیت. ساختمان‌ها و شریان‌های ارتباطی به طور ویژه آسیب پذیرتر می‌باشند. این پدیده 
می‌تواند به خیابانه؛ پل‌ها و بزرگراه‌ها سیب زده حطوط آبرسانی؛ گاز و فاضلاب را مختل کرده به پی ساختمان‌ها 
آسیب رسانده و موجب ترک در آنها گردد. در این میان, سازه‌هایی که وسعت زیادتر و ارتفاع بیشتری دارند سیب 
پذیرترند. بعنوان مثال. حطوط راه آهن. سدهای خاکی. تصفیه خانه‌ها و کانال‌ها از آسیب پذیری زیادتری برخوردار 
هستند. شناسایی مناطق فرونشسست و اندازه‌گیری میزان آن می‌تواند نقش موثری در مدیریت این پدیده و مدل سازی 
آن به منظور پیش بینی ایفا نماید (تورانی و همکاران ۱۳۹۷). با وجود انجام مطالعات فراوان درباره فروند ت در 
محدوده دشت‌های ایران (جنت و همکاران ۱۳۸۸؛ شریفی‌کیا ۱۳۹۱؛ افضلی و همکاران, ۱۳۹۲؛ صالحی و همکاران, 


۲ یر شتاهی رن هیکارا۵ :#۱۲۹۲ بای مرش ۵ ۱۳اب همکارانم ۱۳۹۷ فروشست بهتاهای 
شهری کمتر مورد توجه قرار گرفته است بطوری که می‌توان تنها به مطالعات محدودی اشاره کرد به طور مثال کریمی 
و همکاران (۱۳۹۲): خطرپذیری سکونتگاه‌های شهری واقع در منطقه ۱۸ شهر تهران ناشی از پدیده فرونشست مورد 
بررسی قرار داده و نرخ فرونشست را تا ۲۵ سانتیمتر در بخش‌های جنوبی محدوده مورد مطالعه خود برآورد می‌کنند. 
توزاتن ورمکا ۱۳۹۷0 با استقاوه از روش ناخ سس رادار یه هنت فرو تست ور قرب اسان کلستان 


بویژه شهر گرگان پرداخته‌اند» ایشان محدوده فرونشستی به میزان ۶/۸ سانتیمتر در شمال شهر گرگان مشخص کردند. 


مصص 1 


۱۷۰ جغرافیا و مخاطرات محیطی شماره سی‌ام 


لیو و همکاران! (۲۰۱۷) با استفاده از تکنیک پراکنش گرهای دائمی تصاویر رادار تغییرات ارتفاعی سسطح زمین در 


شهر کانچو چین را مورد بررسی قرار داده و مناطقی که دچار فرونشست و يا بالاآمدگی شده‌اند را مشخص کرده 


است. کاستلازی و همکاران" (۲۰۱) با استفاده از روش تداخل سنجی بر روی تصاویر رادارست» نرخ فرونشست 
شهر مکزیکوسیتی را بطور متوسط سالانه ۱۰ سانتیمتر برآورد کردند. طبق نتایج ایشان مناطق مرکزی و شرقی شهر 
بیشتر دجار فرونشست شده است. 

شهر کرمان نیز به مانند بسیاری از شهرهای کویری ایران به دلیل خشکسالی‌های چند سال اخیر مجبور به استفاده 


بی‌رویه از منابع آب زیرزمینی شده است. در نتیجه با خطر فرونشست روبرو است. در این رابطه عباس نژاد و حسن 


زاده (۱۳۸۵) با استفاده از عکس‌های هوایی و اطلاعات نقشبرداری به ارزیابی خحطر فروندٌ ِ و نورم زمین در 
سطح شهر پرداخته‌اند نتایج مطالعات ایشان نشان می‌دهد که در محدوده‌هایی از سطح شهر فرونشست اتفاق افتاده 
است همچنین ایشان میزان بالاآمدگی را در بخش‌هایی از شهر بالغ بر 1 سانتیمتر اندازه گیری کرده‌اند. ولیکن در 


تحقیق حاضر سعی شده تا با استفاده از روش نوین تداحل‌سنجی تصاویر راداری ۸۵٩۸‏ و 1 ,5۳1101۳1 مقوله 


3 روند 


فرونث زمین در محدوده ‏ کرمان مورد بررسی جامع قرار گرفته. و ضمن تعیین میزان فرو 
تغییرات آن در شهر کرمان طی چهار دوره و در مجموع در طول ۱۵ سال اخیر پایش شده و در نهایت محدوده‌های 


۲- مواد و روش‌ها 

۱-۲- منطقه مورد مطالعه 

شهرستان کرمان با وسعتی در حدود 110۰۰ کیلومتر مربع در بخش شمال شرقی استان کرمان قرار گرفته است. 
این شهرستان از شمال به خراسان جنوبی, از شرق به بیابان لوت. از غرب به شهرستان رفسنجان و از جنوب به 
شهرستان بم و بردسیر منتهی می‌گردد. شهر کرمان در بخش غربی شهرستان و مابین طول جغرافیایی ۵۷ درجه تا ۵۷ 
درجه و ۷ دقيقه شسرقی و عرض جغرافیایی ۳۰ درجه و ۱۶ دقيقه تا ۳۰ درجه و ۱۹ دقیقه شسمالی واقم شده است 
(شکل .)٩‏ ارتفاع متوسط شهر کرمان از مسطح دریا ۱۷۳ متر بوده و دارای آب و هوای خشک و نسیتاً گرم است 
(عباس نژاد. ۱۳۸۳). 


1 10 ۰ 
2 2510112271 6 ۰ 


شکل ۱- موقعیت محد وده مورد مطالعه 


در این پژوهش از دو سری داده راداری استفاده شد. تعداد ٩‏ تصویر از سنجنده ۵5۸۴ ماهواره 0۲1586 سازمان 
فضایی اروپا با فرمت ۹17 (تک منظر) از نوع مد 10886 با پلاریزاسیون ۷۷ که مربوط به بازه زمانی ۲۰۰۶/۰۳/۰۹ 
تا ۲۰۱۲/۰۳/۲۷ است؟؛ و ۲ تصویر از سنجنده 1 ,617۳11۳1 ماهواره 5072 سازمان فضایی اروپا با فرمت 516 از 
نوع مد 10896 با پلاریزاسیون ۷ مربوط به تاریخ‌های ۲۰۱6/۰۷/۰۲ و ۲۰۱۷/۰۹/۲۱. این نوع داده‌ها داده‌های خامی 
هستند که تنها پردازش اولیه بر روی آنها انجام گرفته و به تصاویر تک منظر تبدیل شده‌اند. و اطلاعات آنها به هیچ 
عنوان مخدوش نگردیده است. ترکیب داده‌های راداری از مدارهای صعودی و نزولی " نیز جهت بهبود مدل‌های 
رقومی زمین یا به دست آوردن جابجایی‌ها در جهات مختلف با استفاده از اینترفرومتری می‌تواند مورد اسستفاده قرار 
گیرد. با توجه به اينکه در این تحقیق هدف پایش میزان فرونشست است. تنها از داده‌های نزولی استفاده شده است. 
همچنین مدل ارتفاع رقومی ۳۰ متری 518171 برای منطقه مورد مطالعه استفاده شد. جدول شسماره ۱ مشخصات 


داده‌های مورد استفاده در این پژوهش را نشان می‌دهد. 


حدول ۱- مشخصات داده‌های مورد استفاده در تحقیق 


تست نیع تصویر تاریخ فریم مسیر | خط مبنا عمودی (متر) 
۱ ۰۳-4 ۹ ۲۰۰۷/۰/۲۱ ۳۹۹۷ ۱3۳ ۱۵ 
هد ۱ 2۲ ۲۰۱۰/۰۹/۲۸ ۲۳۹۹۷ ۱۳ 1 
۳ ۳-۹ ۷ 2۲ ۲۰۱۲/۸/۲۶ ۳۹۹۷ ۱1۳ ۱۸ 
3 1 ,5۳۲1۳۲ ۷۲ ۲ ۲۰۱۷/۰۹/۲۱ ۳۰۹۹ ۳۰۹ ۲۷ 
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ام جغرافیا و مخاطرات محیطی شماره سی‌ام 


تحقیق حاضر با استفاده از روش تداحل سنجی راداری صورت گرفته است. این روش در میان روش‌های زمینی 
و فضایی به‌عنوان کارآمدترین روش برای اندازهگیری تغییرات سطح زمین با دقت و قدرت تفکیک مکانی بسیار با 
به شمار می‌رود (2007 .21 6 تلاعتهل. در تداخحل سنجی راداری؛ فاز تصویر ۳ شده از موقعیت‌های تصویر 
تفاضل‌گیری بین این مقادیر تصویر جدیدی حاصل می‌شود که تداخل‌سنج" نام دارد. تشکیل تداخل سنج یا 
اینترفرو گرام اساس پردازش راداری اس (2012 ,۳1006). بدین منظور با استفاده از ترم افزار 60 و ۲۸۳۲( 
به پردازش تصاویر دو زمانه راداری مربوط به سنجنده ۱۳۵ و 1 ٩121110۳‏ با روش فوق اقدام گردید و 
حاصل این فرآیند با بهره گیری از محیط نرم افزاری 103 015 ۸۲۰ علاوه بر اسستخراج نرخ و دامنه فرورفتگی و 
بالاآمدگی, پایگاه داده‌های مکانی ایجاد و پهنه‌بندی این پدیده در منطقه مورد مطالعه صورت گرفت. 


۲-۲- تداخل سنجی راداری 

نحوه اندازهگیری تغییرات سطح زمین به روش تداخل سنجی راداری در شکل شماره ۲ نمایش داده شده است. 
در این شکل 7 معرف فضایی مشخص در سطح است که در قالب یک پیکسل تصویر شده است. سنجنده تصویر 
نحست آ این فضا را در زمان ما ثبت و مقدار فاز آن را اندازه گیری می‌کند (1۷10). مقدار نشست فاصله ( تا ۳1 
است؛ که طی زمان مشخصی صورت گرفته است (106). برای اندازه گیری این مقدار سنجنده تصویری دومی " در 
زمان ) و با هندسه‌ای کاملاً سبیه به تصویر نخست اخذ کرده مقدار فاز را بر روی آن اندازه گیری می‌کند (5)). 
روش تداخل سنجی تفاضلی تفاضل فاز 0 و 1۷1 را در فرم تداخل نگار فازی نمایش می‌دهد (2۳1۳6). در 
صورت ابت و پایدار بودن سطح. تفاضل فاز این دو تصویر (1۷1۳ -[5) ناشی از تغییر موقعیت این دو سنجنده است 
و مقدار آن به وسیله رابطه زیر بدست می‌آید: 

۱ ره - #«نه۵ 

در اندازه گیری مقدار بالاآمدگی به کمک این روش. سطح ناپایدار تصور شده به طوریکه سطح از ۴ به ۳1 تنزل 
یافته است. تعیین مقدار جابجایی در فاصله زمانی دو تصویر (0) تابع احتلاف فاز دو تصویر (۸۵06) به همراه فاز 
ناشی از اثر توپوگرافی (2002)) و فاز ناشی از جابجایی سطح (0[۷۲0۷)) و همچنین فاز ناشی از اثر اتمسفر (۲۸) 


7 
۵:۲ < 7 7 


/ 
< 4100 + ۵۷0۷ + ۰ ۲ 


( 0 
2 96۲ 
3 6 


سال هشتم پایش روند تغییرات ارتفاعی سطح زمین در شهرکرمان و ... ۱۷۳ 


در این روش در صورت در اختیار نداشتن تصویر سوم به کمک مدل رقومی ارتفاعی زمین و تبدیل ارتفاع به 
فان یک تداخل نگار مصسنوعی تولید می‌نسود و از این راه به کمک معکوس اطلاعات 19۳۷ اثر فاز نانسی از 
توپوگرافی محاسبه و از مقادیر احتلاف فاز حذف می‌شود. اختلاف فاز باقی‌مانده به اثر جابجایی سطح و اتمسفر تعلق 
دارد. در نهایت با نادیده انگاشتن اثر اتمسفر در جابجایی‌های به میزان بالا (چند سانتیمتر) و يا حذف آن به کمک 
تصاویر اپتیکی. اختلاف فاز دو تداخل نگار (۸(101/ظ)) (تداخل نگار اصلی و مصنوعی) فقط بیان کننده مقادیر 
جابجایی سطح خواهد بود (شریفی کیاء ۱۳۹۱). فلوچارت روش تداخل‌سنجی راداری در شکل شماره ۳ ارائه شده 


اسیت: 


شکل ۲- الف - هندسه روش تداخل سنجی. ب- روش تداخل سنجی تفاضلی برای اندازه گر تغییرات سطح 
زمین (2000 ,21 )6 مصحصع ۳۳7 


۱ تهیه تصاویر ۸۹۸۴ أ 


۰ 


ماس بسن لاین و انتخامی 
تصاویر مناسي 


ّ 


أ لهیه تداخل لگاشت أ 


۱ حذف فاز ثوپو گرافی با استفاده از ۱ 


سدل ارتفاعی رقوسی <1212۳) 


1 


۱ خاهش نویر یا فیلت رگذاری ۱ 


/ 


۱ واپمچشی فاز (۱ وه ل) ۱ 


شکل ۳- فلوجارت روش تداخل سنجی راداری 


با تهیه زوج تصاویر مناسب (جدول ۱) یعنی زوج تصاویری که خط مبنای عمودی آنها بسیار کم است. اقدام به 


اجرای روش تداحل‌سنجی راداری با استفاده از نرم‌افزار 5 0( در محرط 5 شد. بدین ترتیب پس از 


۱۷ جغرافیا و مخاطرات محیطی شماره سی‌ام 


انجام فرایند بازیابی فازن به منظور رفع خحطای نوفه با انتخاب نقاط کنترل زمینی (660۳) فاز ویرایش شسد. سپس به 
منظور تفکیک سیگنال تغییر شسکل. مژلفه توپوگرافی با استفاده از مدل ارتفاع رقومی» تصحیح و پالایش و تداخحل 
نگاشت‌های مربوطه تهیه شد (شکل 4). به منظور شناخت محدوده‌های پرخطر در سطح شهر کرمان اقدام به تهیه 
نقشه الگوی فضایی عرصه‌های پر مخاطره گردید بدین منظور دو فاکتور نرخ و دامنه فرونشست و بالاآمدگی 
مستخرج از تداخل نگاشت مربوط به سال‌های ۱۳۹۳ تا ۱۳۹۲ به‌عنوان فاکتور عامل در زایش مخاطره به تحلیل گر 
سیستم در محیط 015 وارد شد. نقشه خروجی معرف سلول‌های واجد ارزش متفاوت از فاکتور زایش مخاطره است. 
بر این اساس دسته‌بندی ارزش سلول‌های معرفی شده در نقشه خروجی به ۶ کلاس متفاوت برای خطر فرونشست و 
۳ کلاس متفاوت برای خطر بالاآمدگی. تحت عنوان پهنه‌های برخوردار از مخاطره فرونشست و بالاآمدگی مشتمل بر 
پهنه بسیار پرخطر پهنه پرخطر» پهنه نسبتاً پرخطر و پهنه نسبتاً کم خطر تهیه گردید. 
۳- نتایج و بحث 

همان‌طور که در شکل شماره ۶ مشاهده می‌شود در تصاویر تداحل نگاشت بدست آمده فرینج هایی تشکیل شده 
است. با توجه به نحوه قرار گیری سیکل‌های رنگی مکانیزم جابجایی نیز متفاوت است. بطوری که اگر سیکل رنگی 
زرد-آبی-قرمز باشد. جابجایی دور از رادار و اگر سیکل رنگی زرد-قرمز-آبی باشد جابجایی به سمت رادار اتفاق 
فتاده است. آنچنان که در شکل شماره 4 به خوبی مشسخص است موقعیت فرینچ‌های ایجاد شده دقیقاً از مرز شهر 


تبعیت کرده که نشانگر تغییرات ارتفاعی متفاوت در سطح شهر است. همچنین روند و شکل سیکل‌های ایجاد شده 
در تمامی کوزه‌ها تسا متنابه اتبیتتا: 


سس ت ۳ ت- ۳ ود 2-۹ 


شکل ۶- تداخل‌نگاشت‌های تشکیل شده ۸: مربوط به سال‌های ۱۳۸۳-۱۳۸۲ و 8: مربوط به سال‌های ۱۳۸۹- 


۲ : مربوط به سال‌های ۱۳۸۹-۱۳۹۱ و 8: مربوط به سال‌های ۱۳۹۳-۱۳۹۲ 


به منظور تهیه نقشه جابجایی سسطحی منطقه باید با اعمال فیلتر «اع0010) بر روی تداخل‌نگاشت‌هاء خطاهای 
احتمالی ازجمله خطای اسپکل, را رفع نمود. از آنجا که اطلاعات فاز یک تداخل نگاشت در مقیاس۲7 اسست یک 
مشکل مبهم در محاسبه تعداد صحیح سیکل‌های فاز وجود دارد که این حل ابهام به‌عنوان حل ابهام فاز ( ۳2256 
0۵۵۵ شناخته می‌شود. به فرآیند بازیابی مقادیر فاز نامبهم و صحبح از یک مجموعه دوبعدی مقادیر فاز که 
در بازه (-7 ,7) است. بازیابی فاز دوبعدی گفته می‌شود؛ و در نهایت با انجام و اعمال فیلتر و حل ابهام فاز, نقشه‌های 
جابجایی رخداده در طول دوره‌های مختلف و الگوی توزیع فضایی آن تهیه شد (اشکال ۵ و ). از آنجایی که 


جابجایی‌های رخداده در راستای دید ماهواره اتفاق افتاده, می‌توان مقادیر مثبت را به معنی نزدیک شدن سطح زمین به 


ماهواره و بالاآمدگی و مقادیر منفی به معنی فاصله گرفتن از ماهواره و فروند ت تلقی نمود. در این صورت 
محدوده‌های قرمز رنگ مناطق فرون:ٌ تَ و محدوده‌های آبی رنگ مناطق بالاآمده سطح زمین را نشان می‌دهد. 
همان‌طور که در اشکال شماره ۵ و 7 مشاهده می‌شود میزان فرونشست و بالاآمدگی سطح زمین بین سال‌های ۱۳۸۳ 


ت| ۱۳۸۲ بین ۱۲-۲ ۱۸ سانتیمتر فروندٌ ِ و بالاآمدگی اتفاق افتاده است. محدوده‌های فرونشستی بیشتر در اطراف 
محدوده شهر کرمان متمرکز هستند درحالی که بالاآمدگی سطحی بیشتر در مناطق مرکزی شهر اتفاق افتاده است. میزان 


بای من ستان‌های ۱۳۸۱ ۱ ۱۱۸۹ کی بعش شم انیت وین وی سا ۵ شاف ورشته است رن اسان 


شکل شماره 7 میزان جابجایی طی دو سال مابین سال‌های ۱۳۸۹ تا ۱۳۹۱ بین ۱۸- تا ۱۰ سانتیمتر بوده است» 
همچنین در طول سال‌های ۱۳۹۳ تا ۱۳۹7 اين تغییرات بین ۳۲- تا ۱۶ سانتیمتر رسیده است. آنچه که در دوره‌های 
مختلف مشاهده می‌شود روند افزایشی میزان فرونشست و روند کاهشی میزان بالاآمدگی در طول این سال‌ها است؛ 
به‌طوری که نرخ فرونشست طی دوره ۸۱-۸۳ از ۷/۳ به ۷/۲ سانتیمتر در سال طی دوره ۸٩-۸۲‏ رسیله است و در 


همین مدت میزان بالاآمدگی نیز از 7 به 7/7 سانتیمتر در سال افزايش یافته است. همچنین طی سال‌های ٩۱-۸۹‏ و 
٩۱-۳‏ نرخ فرونشست و بالاآمدگی به ترتیب از ٩‏ به ۱۰/۹ سانتیمتر در سال افزایش و از ۵ به ۶/0 سانتیمتر در سال 
کاهش يافته است. همچنین همان‌طور که در اشکال شماره ۵ و 7 مشخص است. الگوی فضایی جابجایی‌ها نیز طی 
دوره مطالعه تغییر کرده است. به‌طوری که محدوده‌های فرونشست افزايش وسعت نشان داده و تقریباً از اطراف شهر 
کرمان به درون محدوده‌های شهری منتقل شده است در حالی که وسعت مناطقی که دچار بالاآمدگی شدهاند. ابتدا 


یک روند افزایشی داشته ولی از سال ۱۳۸۹ به بعد روند کاهشی در پیش گرفته‌اند. 


۱۳ 


شماره سی‌ام 


3 وم سیر مه ت- سیم امس یه 


شکل ۵- نقشه‌های جابجایی ۸: مربوط به سال‌های ۱۳۸۲-۱۳۸۳ و 8: مربوط به سال‌های ۱۳۸۹-۱۳۸۲ 


شکل *- نقشه‌های جابجایی ۸: مربوط به سال‌های ۱۳۹۱-۱۳۸۹ و 8: مربوط به سال‌های ۱۳۹۲-۱۳۹۳ 


شکل شماره ۷ محدوده‌های پرخطر در سطح شهر کرمان را نمایش می‌دهد. استخراج عرصه متعلق به اين پهنه‌ها 
نشان می‌دهد که قریب ۷۰۰ هکتار از محدوده شسهرکرمان در پهنه‌های پرخطر و بسیار پرخطر به لحاظ خطر 
فرونشست و قریب به ۱۷۰۰ هکتار از محدوده شسهر در پهنه‌های پرخطر و بسیار پرخطر به لحاظ خطر بالاآمدگی 
سطح زمین» واقع شده است (جداول ۲ و ۲). 


جدول ۲- مساحت و درصد عرصه‌های شهری واقع در پهنه خطر فرونشست 


عرصه خطر مساحت (هکتار) درصد 
بسیار پر خطر 1 ۳۸ 
پرخطر ۳ ۳/۸۱ 
نسبتاً پر خطر ۱1۰۹ ۸۷۵ 
نسبتاً کم حطر 9۸۸۰ ۳/۹۷ 


جدول ۳- مساحت و درصد عرصه‌های شهری واقع در پهنه خطر بالاآمدگی 


عرصه خطر مساحت (هکتار) درصد 
تیاه پزشتظر ۳۱ ۳/۷ 

پرحطر ۱۳۱ 1/۳ 
نسبتاً پرخطر ۲۳۷۸ ۱۳/۲ 


سال هشتم پایش روند تغییرات ارتفاعی سطح زمین در شهرکرمان و ... ۱/۷ 


گِ- 
۳ 


شکل ۷- پهنه‌بندی مخاطره فرونشست (۸) و بالاآمدگی (8) سطح زمین در مناطق چهارگانه شهر کرمان 


با روی‌هم‌اندازی لایه پهنه‌های پرمخاطره بر روی نقشه مناطق چهارگانه شهرکرمان. محدوده‌های شهری 
پرمخاطره دقيقاًقابل شسناسایی می‌شود بطوری که طبق شکل شماره ۷ سه محدوده از شهر کرمان شامل محدوده 
بلوار پارادیس, منطقه حول میدان مشتاق در بخش مرکزی و بافت تاریخی شهر و شهرک صنعتی در بخش شرقی 
شهر جزء مناطق پرمخاطره به لحاظ فرونشست قلمداد می‌شوند. همچنین به لحاظ مخاطره بالاآمدگی سطح زمین نیز 
۲ محدوده جزء مناطق پرمخاطره محسوب می‌شوند. اين مناطق شامل محدوده بلوار سیدی در بخش شمالی شهر و 
شهرک مطهری در بخش جنوبی شهر می‌شود. مطالعات میدانی انجام گرفته در تمامی محله‌های ذکر شده نشانگر 
تأثیر زیاد این مخاطرات بر مناطق مسکونی می‌باشد. به‌طوریکه شواهد زیادی مانند ایجاد ترک در دیوارهای بعضی 
از منازل» نشست سنگفرش‌ها و تخریب تأسیسات شهری قابل مشاهده است (اشکال ۸ و .4٩‏ 


شکل ۸- نمونه‌هایی از تاثیرات فرونشست بر سازه‌های شهری 


اش جغرافیا و مخاطرات محیطی شماره سی‌ام 


شکل -٩‏ نمونه‌هایی از تأثیر بالاآمدگی زمین بر ساختمان‌های شهری 


نتایج به دست آمده نشان می‌دهد که روند کلی تغییرات سطح زمین در محدوده شهر و اطراف آن دو روند کاملاً 
معکوس بوده است. بدین ترتیب که مناطق پیرامونی شهر در حال فرونشست بوده ولی در محدوده‌های شهری سطح 
زمین دچار بالاآمدگی شده است و این روند با کمی تغیبرات هم‌اکنون نیز ادامه دارد با این تفاوت که در چند سال 
اخیر بخشی از مناطق شهری نیز دچار فرونشست شده‌اند و روند بالاآمدگی سطح زمین نیز بسیار کاهش یافته و تنها 
به بخش‌های کمی محدود شده است. دلیل تغییرات رخ‌داده در ارتفاع سطح زمین در طول این چند سال را می‌توان 
این گونه تفسیر کرد که افزایش آب شرب مصرفی که ناشی از افزايش جمعیت و توسعه شهر کرمان می‌باشد. همراه با 
انسداد مجاری زیرزمینی مانند قنوات. تغییر کاربری اراضی کشاورزی» حذف جاه‌های بهره برداری و استفاده از 
به‌طوری‌که این استدلال با نتایج مطالعات آقاملایی و همکاران (۱۳۹۲)» نیز همخوانی دارد. علاوه بر آن به دلیل نوع 
خحاک شهر کرمان که بیشتر از نوع خاک‌های رسی می‌باشد (مینی‌زاده و همکارا ۱۳۹۵) در نتیجه جذب آب توسط 
خاک‌های غنی از رس باعث تورم و بالآمدگی سطح زمین شده است. این مسئله در نوع خرابی ساختمان‌ها نیز 
مشهود است همانطور که در شکل شماره ٩‏ مشاهده می‌شود. پی بسیاری از منازل واقع در محدوده‌های بالاآمده دجار 
خیس خورد گی ناشی از بالا بودن سطح آب‌های زیرزمینی شده است. البته با احداث سیستم فضلاب در چند سال 
اخیر روند بالآمدگی سطح زمین نیز کاهشی شده و این موضوع در نتایج کاملاً مشهود است؛ و اما حفر چاه و برداشت 
شهر کرمان شده است؛ و دلیل گسترش این مناطق فرونشستی به محدوده شهر کرمان می‌تواند علاوه بر حفر چاه 
مقایسه سطح برخورد به آب زیرزمینی (۱۳۹۵-۱۳۵) در دشت کرمان و محدوده شهر کرمان توسسط آقاملایی و 
همکاران (۱۳۹۲) نشانگر معکوس بودن این دو روند است. به گونه‌ای که در برابر افت سالانه سطح آب در دشت. با 


افزایش سالانه سطح آب در محدوده شهر مواجه بوده‌ايم که نتیجه آن در تغییرات ارتفاعی سطحی رخداده کاملاً 


سال هشتم پایش روند تغییرات ارتفاعی سطح زمین در شهرکرمان و ... ۷۹ 


مشهود است. مقایسه نتایج به دست آمده از این تحقیق با مطالعات عباس‌زاده و حسن‌زاده (۱۳۸۵), نشان می‌دهد که 
مناطق فرونشستی روند تغییرات مشابهی را طی نموده است. علاوه بر آن ایشان. میزان بالاآمدگی زیاد سطح زمین 
اندازه‌گیری شده توسط دستگاه‌های نقشه‌برداری را ناشی از حطای احتمالی این تجهیزات دانسته است در حالی که 
نتایج این تحقیق آن را واقعی دانسته و تأیید می‌کند. مقایسه نتایج مطالعات مشابه ( 64 هنتته[[عاهه ,2017 .2 1 دازآ 
6 .۵) با نتایج حاصل از این پژوهش نشانگر یکسان بودن دلایل زمینه‌ها و عوامل رخداد و وقوع اینگونه تغییرات 
در سطح زمین. که همانا تغییرات سطح آبهای زیرزمینی بوده است. 


۶- جمع بندی 

امروزه تصاویر راداری در بسیاری از مطالعات جغرافیایی کاربرد دارد. یکی از موارد کاربرد آن در تعیین 
جابجایی‌های ارتفاعی سطح زمین است امری که در این تحقیق از آن استفاده شده است. نتایج بدست آمده از روش 
تداحل‌سنجی راداری نشان می‌دهد که مناطق مختلف شهرکرمان در طول ۱۵ سال اخیر به لحاظ فرونشست و 
بالآمدگی رفتارهای متفاوتی داشته به‌طوری که در طول سال‌های ۱۳۸۲-۱۳۸۳ با توجه به تمرکز اراضی کشاورزی و 


به تبع آن افت نسبی سطح آب‌های زیرزمینی در اطراف و خارج از محدوده شهر, مناطق فرونشستی بیشتر در اطراف 


شهرکرمان متمرکز بوده و نرخ فرونشست نسبتاً ایین بوده است. همچنین در اين دوره به دلیل باآمدگی سطح 
آب‌های زیرزمینی در محدوده شسهر, بخش وسیعی از مناطق مرکزی شهر دچار پدیده بلآمدگی سطح زمین شده 
است. بین سال‌های ۱۳۸۹-۱۳۸۲ موقعیت مناطق فرونشست در اطراف شهر کمی تغییر کرده و محدوده‌هایی از شهر 
نیز درگیر پدیده فرونشست شده است. در طول سال‌های بعد مناطق بیشتری از سطح شهر دچار فرونشست شده و 
تاکنون نیز ادامه دارد. پایش بالاآمدگی سطح زمین طی این ۱۵ سال کمی متفاوت بوده است به‌طوری که طی سال‌های 
۳ تا ۱۳۹۲ روند گسترش مناطق بالاآمده کاهشی شده که این می‌تواند ناشی از پایین رفتن سطح آب به‌دلیل 


راه‌ندازی و احداث سیستم فاظلاب شهری و در نتیجه پایین آمدن رطوبت خاک‌های رسی و نشست زمین باشد. با 


توجه به وجود شواهد تخریبی در سازه‌های موجود در محدوده‌های بالاآمده هر می‌توان چنین برداشت کرد که 
باآمدگی سطح زمین نیز به مانند فرونشست می‌تواند مخاطره آمیز باشد. با توجه به پهنه‌بندی حطر فرونشست شهر 


کرمان. سه منطقه شامل محدوده بلوار پارادیس منطقه حول میدان مشتاق در بخش مرکزی و بافت تاریخی شهر و 
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